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The foundational theory for dusty plasmas is the dust charging theory that provides the dust potential and 
charge arising from the dust interaction with a plasma. The most widely used dust charging theory for 
negatively charged dust particles is the so-called orbital motion limited (OML) theory, which predicts the 
dust potential and heat collection accurately for a variety of applications, but was previously found to be 
incapable of evaluating the dust charge and plasma response in any situation. Here we report a revised OML 
formulation that is able to predict the plasma response and hence the dust charge. Numerical solutions of 
the new OML model show that the widely-used Whipple approximation of dust charge-potential relationship 
agrees with OML theory in the limit of small dust radius compared with plasma Debye length, but incurs 
large (order-unity) deviation from the OML prediction when the dust size becomes comparable with or larger 
than plasma Debye length. This latter case is expected for the important application of dust particles in a 
tokamak plasma. 


I. INTRODUCTION 

Two years before Irving Langmuir coined the term 
“plasma” as in plasma physics^, he and Mott-Smith laid 
down a foundational theory on the charging of a spheri¬ 
cal and cylindrical probe in a laboratory plasmai^, which 
was necessary for interpreting the measurement of what 
became later known as Langmuir probes^. In the ensuing 
decades, the charging and dynamics of solid particulates 
immersed in plasmas built the foundation for a new dis¬ 
cipline in plasma physics - dusty plasmas^^— , in which 
the collective behavior of a group of dust particles in 
a plasma environment is studied. The physics of both 
dust in a plasma and dusty plasmas finds applications in 
space^i^, astrophysics^iiSj and laboratoryii^—. The fun¬ 
damental dust-plasma interaction, which is also essential 
to understand the collective behavior in a dusty plasma, 
includes (1) charging of the dust by absorption of plasma 
particles so it becomes subject to the electric field (E) via 
the electrical force Q^E {Qd is the dust charge); (2) heat¬ 
ing of the dust via the collection of plasma electron and 
ion energy fluxes so the dust particulate can melt, evap¬ 
orate, or simply sublimate; (3) dragging of the dust via 
a frictional force enhanced by the Coulomb interaction 
with the flowing background plasma. 

To understand the dust dynamics and its change of 
state, one must resort to a dust charging theory. The 
most-widely used of such is the Orbital-Motion-Limited 
(OML) theory, which dates back to the work of Mott- 
Smith and Langmuir—. In the 1960s, Al’pert et ali^ and 
Laframboisei^ completed the current formulation. The 
OML theory is known to predict accurately the dust po¬ 
tential for a smalli^i2£ and not so small dust grailiii"— , 
despite its simplifying assumption on collisionless ion or¬ 
bit that misses the absorption radius effect away from 
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the dust surfac o ^i^. Surprisingly, the existing formula¬ 
tion can not be used for predicting the plasma response, 
namely solving the plasma potential </>, as shown in Allen 
et al2^. As the result, OML theory to date can not pre¬ 
dict the fundamental quantity of dust charge. In practice, 
one has been using an idealized dust charge-potential 
relationship due to Whipple^, who computed the dust 
capacitance using the conventional Debye shielding po¬ 
tential that does not take into account the constraint of 
angular momentum conservation in setting the plasma 
density near the dust surface. 

The purpose of this paper is to present a revised OML 
formulation that, for the first time, is able to predict 
the plasma response and hence the dust charge. The 
OML predictions will then be contrasted with the Whip¬ 
ple approximation to elucidate the missing physics in 
dust charging for a spherical dust of comparable size with 
the Debye length, which is of special importance to dust 
transport/survivability in tokamak a^^’^^ . This new OML 
formulation is also important for electron-emitting dust 
particulates as it provides the basis for extending the 
OML theory to positively charged dusti^^. 

As an approximation of the Orbital Motion (OM) the¬ 
ory which follows the collisionless particle orbit via con¬ 
servation of energy and angular momentum, the OML 
theory will inevitably introduce discrepancies in the eval¬ 
uation of dust current/heat collection, and in the plasma 
ion density evaluation. The former will give rise to a dis¬ 
crepancy in dust potential, while the latter in plasma po¬ 
tential and hence dust charge, in comparison with those 
predicted by the OM theory. The usefulness of the OML 
formulation is due to its simplicity, and the relatively 
high accuracy for a variety of applications where the dust 
size is not large compared with the plasma Debye length. 
It must be clarified that the fatal breakdown of the pre¬ 
vious OML formulation in the form of an imaginary ion 
density when applied to calculate the plasma potential 
and dust potential^^ is not due to the OML approxima¬ 
tion of ignoring the ion absorption radius effect. The root 
cause is the OML ion density formula originally given in 



Ref which we find will produce an imaginary density 
as long as Debye shielding of a dust particle is in place. 
In the appendix, after we present the corrected OML 
theory and its application in calculating the plasma po¬ 
tential and dust charge, we give a detailed account of how 
the previous OML formulation would break even in the 
regime of small dust size compared to Debye length, a 
limit commonly accepted for OML approximation to be 
highly accurate. 

The rest of the paper is organized as follows. In sec¬ 
tion El we briefly recall the key components of the OML 
theory and its inability to predict dust charge in its cur¬ 
rent formulation. A corrected formulation on the OML 
ion density is shown in section uni which completes the 
OML formulation for calculating both dust potential and 
dust change. In section llVl we apply the new OML model 
to evaluate the plasma response and the dust charge, and 
contrast the results with the widely used Whipple ap¬ 
proximation. The findings are summarized in section [V] 



FIG. 1. The effective potential <I>efi is shown for ions 
in a hydrogen plasma with different values of J or a = 
/{2mikBTe). $eff can have extrema at Vm > rd that re¬ 
flects some ions (a = 2.5, 3.0). The OML approximation ig¬ 
nores this absorption radius effect, which can over-estimate 
the ion current {a = 2.5) and the ion density between dust 
surface and the absorption radius (q = 2.5, 3.0). 


II. BACKGROUND ON OML THEORY 


There are three essential ideas underlying OML the¬ 
ory. First, the collection of plasma electrons and ions by 
the dust is governed by their collisionless orbit. These 
are subject to two conservation laws: energy [E = 
'nia {vr + /2 + Conservation and angular mo¬ 

mentum (J = mvtr) conservation. Here m (g) and (j) are 
the mass (charge) of the particles of species a {a = e, i 
labels electrons and ions, respectively) and the plasma 
potential, while r, Vr and Vt are the radial distance, radial 
velocity and tangential velocity in a spherical reference 
frame centered on the dust grain. The radial motion 
of the plasma particles is governed by a one degree-of- 
freedom Hamiltonian with effective potential $6117 

1 1 

H = -niavl + ^es{r), $eff = T:-J + 9a(A- (1) 

2 2 niar^ 

Second, whether the plasma particle reaches the dust 
is governed by For a typical negatively charged 

dust, 4> is negative and monotonically increasing with 
r, so that $611 is positive and monotonically decreas¬ 
ing with r for electrons. As the result, only electrons 
with E > = I'd) from far away can reach the 

dust of radius r^. The final and third idea, which is a 
simplifying approximation, is what delineates OML from 
the more complete but rarely used Orbital Motion (OM) 
theoryi^i^. The OML theory approximates the ^eff for 
ions also as a monotonic function of r, which results in the 
simplification that ions with > —2^es{'rd)/fni from 
far away can reach and charge the dust. This is a re¬ 
markable simplification since one does not need <f)(r) to 
evaluate the ion and the electron current collected by the 


dust, which are 


le = 

li = 


-eAirrlueo^ 

Ze4:TTrdZnio 


ksTe ( e4>d \ 

27rme yksTe) ’ 


1 ksE 

/ Ze(l)d\ 

V 2'Kmi 

\ keE J 


( 2 ) 

( 3 ) 


Here e is the elementary charge, ks is the Boltzmann con¬ 
stant, rieo {nio) is the electron (ion) density away from 
the dust, rieO = Zriio, is the electron (ion) tem¬ 

perature, qe = —e, Qi = Ze, and 4>d = </>('''d)- Setting 
li + le = 0, one can solve for the dust potential as a 
function of T^/Ti, ion charge state Z, and rrielrrii. For 
an electron-proton plasma with Tg = Ti, OML predicts 
(j)d = —2.5{kBTf,/e), which is in remarkable agreement 
with particle-in-cell simulations^^ that do not make the 
OML assumption of ion being monotonic. 

The OML approximation of a monotonic ion $eff is 
known to be violated for ions with a certain range of J. 
This comes about^i because the centrifugal potential en¬ 
ergy J^/(2TOir^) decreases with r at precisely l/r^, but 
the electrical potential energy Zecj){r) increases at a faster 
rate (exponential in r) in the Debye shielding region, be¬ 
fore it eventually asymptotes to 1 for large r. This can 
produce one or multiple extrema in $efi away from the 
dust surface, Fig.[l] which can turn back some ions in a 
range of J at r™ > r^. In the literature this is known as 
the absorption radius (at > rd) effect for certain ions. 
By neglecting this subtlety, OML approximation would 
over-estimate the ion current to the dust. Interestingly, 
for dust size small and even comparable to Debye length, 
this correction appears to be small and OML prediction 
of dust potential remains reliable. 

A far more serious problem was identified by Allen, An- 
naratone, and de Angelis in 20002^ that the OML theory 
can not predict the plasma response to the presence of 
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a dust particle, which requires the solution of the OML Poisson equation for the plasma potential 

= -£o"^ {Zenf^^ - , (4) 

where Eq is vacuum permittivity. Allen et ali^ used 
the well-known OML ion and electron density given by 
Al’pert, Gurevich, and Pitaevskii in 1965^, 
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where 


z = r/vd, 

(7) 

p = ecji/ksTe, 

(8) 

13 = Te/T,, 

(9) 

<P= {‘P- V’dlz'^) I{1- z-"^). 

(10) 


The failure of OML theory manifests in an imaginary ion 
density when cj) > 4>dlz^. For the initially faster than 1/r^ 
increase in (j){r) due to Debye shielding, which is the cause 
for ion absorption radius at > ra as noted previously 
for a negatively charged dust, Allen et al. concluded 
that the contradiction between (j){r) and an imaginary rii 
in Eq. ([H]) would be an intrinsic defect which prevents 
the OML theory from predicting the plasma response. 
Since the dust charge is related to the normal electric 
field at the dust surface, which requires the solution of 
^(r), one reaches the inevitable position from Ref^ that 
OML theory can not predict the dust charge, despite its 
success on dust potential as re-affirmed in Reb ^^d^ . 


III. RESOLVING THE OML CONTRADICTION 
BETWEEN cf){r) AND ni(r) 

Intuitively it is quite puzzling that the OML simplifica¬ 
tion of ignoring the absorption radius effect, which con¬ 
tributes a small error in the ion charging current, would 
produce an imaginary ion density, as revealed in Allen 
et al.’s analysis. Physically, a maxima in ion $eff(r) at 
I'm > I'd would pose a barrier that turn back ions with 
a range of J. Ignoring this effect with the OML approxi¬ 
mation should manifest in an over-estimation of the ion 
density for r < r^. So why Allen et al. discovered an 
apparently inherent contradiction in the OML theory be¬ 
tween (j){r) and 

Following the physical picture just given, one is 
tempted to conclude that the resolution has to come from 
a revision of the OML expression for as given in 

Eq. dni). We find that this is indeed the case. The cause 
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is a change in integration bound for the OML ion density 
when the plasma potential transitions from ^(r) < (l)dlz"^ 
to (f){r) > (jid/z^. To understand this subtlety, which has 
been evidently elusive for the past five decades, we recall 
that since OML assumes a monotonically varying 4>eff(r), 
the ions at r have a collisionless orbit that will either in¬ 
tercept the dust particle or be reflected by the effective 
potential before it can reach the dust surface. 

In the canonical case that lim^^oo = Oj the birth 
energy of the background plasma ion far away must have 

Eo = {y'^ v'^) -\- Ze(f){r —>• oo) > 0. (II) 

If this ion reaches r, it must have, at r, that 

E = -nil {yl + -I- Ze(f){r) = Eq >0. (12) 

For a negatively charged dust which has cj}{r) < 0, a 
plasma ion of such unbounded orbit (meaning that the 
ion orbit connects to infinity) must have higher kinetic 
energy as it approaches the dust, 

—2Ze4'{r)/mi > 0. (13) 

As illustrated in Fig. [51 this is outside a circle in {vr,vt) 
space, which intercepts the Vr = 0 axis at 

Ut = ^-2Ze<^{r)lrfii. (14) 

Not all of these ions can reach the dust surface (r = r^) 
due to angular momentum conservation. In the OML 
approximation, the effective potential 4>eff('r) peaks at 
r = Td, 

‘^esird) = I + Ze<pd- ( 15 ) 

The ions with E < 4)eff(rd) will be reflected by the ef¬ 
fective potential before they can reach r^. At r > r^, 
these reflected ions satisfy, after explicitly writing out 
E < $eff(rd), 

- {z'^ - l)vt <2Ze{(l)d-(j))/rrii. (16) 
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For a negatively charged dust with (j)d — (t>{r) < 0, the 
reflected ions are bounded by a parabola in Vt) space, 

[z^-l)vl - vl >2Ze{(t)-(fid)/rrii. (17) 


If this is the case, the integration has two zones, sepa¬ 
rated by a set of critical and v^, which are the inter¬ 
cept of the two constraints, v^+vl = —2Ze(f){r)/mi and 
- l) v'^ -vl= 2Ze{(t> - 4>d) /m^, 


As illustrated in Fig. O it intercepts the Vr = 0 axis at 


V 


r 

t 


2Ze 4>- (f>d 
rrii z'^ — 1 


(18) 


As we shall see, this should be compared with the inter¬ 
cept of Eq. (fl^ with tir = 0 axis, i.e. in Eq. da. 


To evaluate the ion density at r > r^, which is the solu¬ 
tion of the Vlasov equation, we integrate the Maxwellian 
distribution, which is assumed for the background plasma 
far away, in the velocity space region where they are al- 
lowedi^. For ions with < 0, their population lies in 
the region of the {vr,vt) velocity space given by 



(23) 


As shown in the left-side diagram of Fig. ^ the inte¬ 
gration can be carried out in two zones, depending on 
whether Vr > or not. The lower zone (denoted as I) is 

Vr € [0,u(:] k Vt G [\/-(2Ze/mi)<p - v^,oo). (24) 

The upper zone (denoted as II) is 


Vr G [u(:, oo) k Vt € 


2Ze 


(0- 




(25) 


Ur < 0 k v1+v1> —2Ze4>{r)lmi. (19) 

For ions with Vr > 0, only the reflected ions are present 
so they are in a region of the velocity space given by 

Ur > 0 k v1 +v1 > —2Ze(j){r)/'mi (20) 

[z^ - l)v1 - vl > 2Ze{(j) - (jid)/rrii (21) 


The ion density is given by 


, , n., ( m. ff 

ni{r) = —= —— / / e vtdvtdvr 

V27r yksTiJ JJi+u 

which yields Eq. (jH), now valid only for <f>{r) < (fd/ 


In the case that u[^ > v\ , 


i.e., 
2 


4> > 4>d/z 


(26) 


These give rise to integration bounds for Ur < 0 and 
Vr > 0 separately. The one for Ur < 0 is straightforward, 
see Fig. [21 and the previous OML formic is correct. 

The one for Ur > 0 is complicated by the possibility 
that Uj ^ < Uj , which implies 

< f>d/z^ (22) 


the integration is restricted to a single zone given by 


Ur G [0, oo) k Vt G 


2Ze 
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which is illustrated in the right-side diagram of Fig. |21 
The ion density then takes the form 


ni{z) 

nto 
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> (j^d/z^, 


(28) 


where (j){z) > (fdjz^, while Eq. ([6|) should be used for 
(/)(z) < (fdjz^. Contrasting Eq. (E51) with Eq. ([51), one sees 
that the previously-known imaginary OML ion density 
for (p > (pd/z'^ is removed in the corrected OML theory 
for Hi- 


IV. OML PREDICTION OF PLASMA RESPONSE AND 
DUST CHARGE 


is important to note that the radial electric held at the 
dust surface Ed = —d(f/dz{z = 1) is related to the dust 
charge according to Gauss’s law, 

Qd = AneorlEd. (29) 

Since the OML theory to date can not predict Qd, users 
of the OML charging theory have been using a simple 
relation for spherical dust capacitance due to Whipple^, 

Qd = 4Treord{l + rd/XD)4>d, (30) 


With the corrected OML ion density, the OML Poisson 
equation ([H) can be solved for the plasma response. It 


with 


Ad = \/ eoTefe^fueo 


( 31 ) 
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FIG. 2. The integration bound in velocity space for the OML ion density, ni{z), changes depending on whether 4i{z) < 
or not. This results in two distinct forms for the ion density in OML theory. 



FIG. 3. The plasma potential (/>, ion density rii, and electron 
density rie are plotted as functions of (r — rd)/\D for three 
cases of Vd/Xo = 1 (black), VdlXo — 0.1 (red), and rd/Xo = 
0.01 (green). 

the electron Debye length. Whipple’s idea is to expand 
the plasma density at z ^ 1, as in the standard Debye 
shielding calculation. Applying this approach to OML’s 
Poisson equation, one finds 

Ue/neO ~ 1 + (32) 

Tii/riioKil- ZI3tp, (33) 


and defining the normalized dust radius as 


d — '^d/ 

(36) 

one finds the solution to the Poisson equation. 


(fi = {(fid/z) exp [-rd{z - 1)]. 

(37) 

Hence the dust charge is 



(38) 

Unlike the original Whipple formula [Eq. (1301)] which as¬ 
sumes a uniform ion density, Eq. (1381) takes into account 
the ion density response. It is interesting to note that for 


a small dust A^, the Debye shielding contribution, 
which is usually small, can be enhanced substantially if 
Z ^ 1 and/or Te^ Ti. We will call the charge-potential 
relationship given in Eq. (l38l) the generalized Whipple ap¬ 
proximation, to distinguish it from the well-known Whip¬ 
ple formula in Eq. (l30l) . 

The actual OML density for r — < Ad can deviate 

significantly from the asymptotic expansion valid when 
r — rd ^ Xd, see Fig. [31 This is already evident from the 
analytical form of rie and in OML theory. Since the 
dust charge can be alternatively computed by integrating 
the net charge of the plasma. 


and hence 


1 d 


dz 


,dip 


dz 






(34) 


Introducing the so-called linearized Debye length 

Xun = Xd/V^ + ZP, (35) 


Qd = - J{Zni-ne)/neod^x, (39) 

due to overall charge conservation, one suspects that the 
OML prediction of dust charge can be significantly dif¬ 
ferent from that in Eq. (1381) . From the numerically eval¬ 
uated Tie and Tii shown in Fig. [31 we find that sharper 
deviation occurs closer to Vd as the dust size becomes 


























smaller, but there is a greater spatial extent of the devia¬ 
tion for large dust size. Since the dust charge corresponds 
to the spatial integration of electron and ion density, one 
finds that such deviation from the simple Debye shielding 
calculation is proportional to rd/Xo- In Fig. |4l we plot 

r = Qd/i'i'!T£ord4>d) - 1 (40) 

as a function of rd/Xo- The specific case has Z = 1 and 
/3 = Te/Ti = 1, so the generalized Whipple approxima¬ 
tion, Eq. (l38l) . is simply 

r = lAUrd/Xo- (41) 

The deviation of the OML prediction of dust charge from 
Whipple approximation can be significantly greater if 
the dust size is approaching the Debye shielding length, 
which is consistent with the results of Refi^. This is 
ultimately due to the effect of angular momentum con¬ 
servation on the plasma density near the dust particle. 
Obviously a significant correction in Qd implies a sub¬ 
stantial change in the electrical force QdE, which can 
impact the dust dynamics, for example, in a tokamak 

reactori22i2i. 


V. SUMMARY 

In conclusion, we have resolved a long-standing issue 
in the OML charging theory, and by doing so, obtained 
a complete OML theory that predicts both dust poten¬ 
tial and dust charge. This is enabled by a revised OML 
ion density formula for the case of </> > (/d/z^i which 
is given in Eq. (E51) . We also provide the first calcu¬ 
lation of the plasma potential and the dust charge us¬ 
ing the OML theory. Our results show that for applica¬ 
tions where the dust particulate radius is much smaller 
than the plasma Debye length, the Whipple approxima¬ 
tion is in good agreement with the OML prediction of 
the dust charge-potential relationship. When the dust 
size is comparable to or larger than the Debye shielding 
length, which is a case of importance to magnetic fu¬ 
sion, we discover significant deviation from the Whipple 
approximation of the dust charge-potential relationship. 
This is attributed to the fundamental role of angular mo¬ 
mentum conservation in setting the plasma electron and 
ion density near the dust particle. 

Since OML is an approximation of the OM theory by 
ignoring the ion absorption radius effect, one expects 
that its prediction of dust potential and charge should 
deviate from that of the OM theory. The usefulness of 
OML charging theory traditionally lies with its surpris¬ 
ingly good agreement with OM on dust potentiaU^^—>2^. 
With the new OML formulation that is able to calculate 
dust charge and plasma reponse, it is an important next 
step to carry out a detailed comparison of the OML dust 
charge prediction with OM theory for a range of dust 
sizes. 



FIG. 4. Dust capacitance is shown as a function of rd/Xo- 
The deviation from Eq. (1381) becomes large when the dust 
size becomes comparable or greater than the Debye length. 
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Appendix A: The breakdown of previous OML forumlation 

As pointed out by Allen et al2^, the OML ion den¬ 
sity can become imaginary if one is to solve the Pois¬ 
son equation, Eq. for the plasma potential using the 
previously known OML ion density formula by Al’pert 
et ali^. This suggests a fundamental breakdown of the 
OML theory for plasma potential and hence dust charge 
calculation. Allen et al investigated this problem from 
the angle of the OML approximation, i.e., the negelect 
of the ion absorption radius effect. In the main text, we 
give a physics argument that ignoring the ion absorption 
radius effect should only introduce a discrepancy in ion 
density compared to OM prediction, but not an unphys¬ 
ical imaginary number. Our corrected calculation of the 
OML ion density removes the possibility of an imaginary 
density. Here we give a more detailed account of how 
the imaginary ion density comes about in the previous 
OML formulation. The objective is to further clarify the 
contrast between (1) a physical approximation (OML) 
that introduces quantitative discrepancy which vanishes 
in its limit of applicability, and (2) an invalid theoretical 
formulation that produces unphysical results. 


The condition for ni{z) to become imaginary (hence 
unphysical) in Eq. ([6]) is to have 


For a negatively charged dust ((/J < 0 and (pd < this 
implies a plasma potential ip{z) that rises faster than 
Pd /with z — 1 the normalized distance from the dust 
surface. As long as Debye shielding is in action (i.e. the 
plasma transport to the dust surface is not intrinsically 
ambipolar), we know that p would rise exponentially as 
a function of 2 within the Debye shielding sphere. The 
cross-over point after which Eq. dAH) is satisfied can be 
illustrated using the Debye shielding potential solution 
of Eq. (l37l) by setting pd /= 1. The cross-over point 
Zc is the solution of 

Zc exp [-fd{zc - 1)] = 1. (A2) 

Here fd is the dust radius normalized against the Debye 
length and z — 1 is the radial distance from the dust sur¬ 
face normalized against the dust size. For = 0.01, one 
finds Zc = 648 or 6.48 Deybe length away. If = 0.0001, 
we have Zc = 116672 or 11.7 Debye length away. In other 
words, the breakdown of the ion density formula, Eq. ®, 
is the direct result of the Debye shielding physics, which 
is independent of how small the dust size is compared to 
the Debye length. 

Of course, the Debye shielding potential as given in 
Eq. (I37II is an analytical approximation, so the actual 
cross-over point Zc would likely differ in its exact loca¬ 
tion from that predicted by Eq. (IA2I1 . When the OML 
Poisson equation is numerically solved, the failure of the 
OML formulation using the ion density given in Eq. ([6]) 
manifests in a not-a-number floating point violation for 
all small r^s we have attempted, as one would expect by 
taking the squared root of a negative number. This is to 
be contrasted with the corrected OML formulation which 
produces the numerical solutions shown in section 


